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We have explored through numerical simulations the amplification of laser pulses in nonlinear fiber Bragg 
gratings with gain and normal dispersion. The distortion of the temporal intensity profile caused by the 
higher-order dispersion coefficients is examined. The stabilizing effect of fourth-order dispersion on pulse 
propagation is demonstrated. For some conditions, the amplified pulse shape evolves towards a Bragg 
similariton, a quasi-parabolic temporal profile with wave-breaking free temporal characteristics. © 2013 
Optical Society of America 
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The amplification of short laser pulses in the normal dis- 
persion regime has been largely studied in recent years 
since it is a key consideration in order to achieve high- 
energy pulses from ytterbium-doped fiber laser systems. 
It is now well known that high intensity pulse propaga- 
tion in an optical fiber with normal group velocity disper- 
sion (GVD) is subject to the optical wave-breaking insta- 
bility. The effects of this instability are known to be alle- 
viated by the use of a parabolic pulse profile [1]. This idea 
has been successfully extended to amplifiers with normal 
GVD by introducing the self-similarity concept [2,3]. It 
has been shown that the nonlinear Schrodinger (NLS) 
equation with gain admits of a self-similar parabolic solu- 
tion (or similariton) characterized by a parabolic tempo- 
ral intensity profile with a linear chirp; the field distribu- 
tion and the pulse spectrum evolve self-similarly during 
amplification, making pulse compression feasible even at 
high power. 

Since similaritons are asymptotic solutions to the NLS 
equation with gain, any arbitrary input pulse will even- 
tually evolve towards the parabolic pulse propagation 
regime. However, in practice, many constraints need to 
be fulfilled by the amplifier and the seed pulses in or- 
der to reach the parabolic regime in a finite length of 
amplifier. It has also been demonstrated that the effec- 
tiveness of the self-similar regime is inherently limited by 
the Raman response of the fiber material and the finite 
bandwidth of the amplifier [4,5]. Therefore, in order to 
obtain the shortest pulses with the largest peak power, 
the length of the amplifier must be kept below a critical 
value. This constraint could be relaxed if one could de- 
sign an amplifier with a larger value of normal GVD. One 
potential avenue to increase the effective GVD would 
consist in propagating the pulses through a fiber Bragg 
grating outside its photonic bandgap. In such a case, one 
could expect to mimic the classical parabolic pulse am- 
plification regime with a much higher value of GVD. 

In this letter, we present the results of a numerical 
investigation of pulse amplification in nonlinear fiber 
Bragg gratings (NLFBG) with gain. The pulses ampli- 



fied in such media, termed Bragg similaritons (BS), have 
been found to share common properties with classical 
similaritons generated in standard optical fibers with 
gain. We will also show that the amplified pulses are 
sensitive to the large value of the third-order dispersion 
(TOD) induced by the fiber Bragg gratings. However, 
our numerical simulations reveal that the fourth-order 
dispersion (FOD) also present alleviates the effects of 
TOD which normally leads to a shock-wave type of in- 
stability. 

The propagation of an optical pulse in an NLFBG with 
gain can be described by the following pair of equations: 
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where Uf and u b represent the field envelopes of 
the electric field distributions corresponding to the for- 
ward and backward propagating waves, v g represents the 
group-velocity in the absence of the grating, 7 is the 
nonlinear coefficient, k,(z) and S(z) are the coupling and 
detuning coefficients of the grating. The constant g rep- 
resents the uniform small-signal gain of the amplifier. In 
an NLFBG without gain or losses, when the peak power 
of the pulse (Pq) is low enough to meet the condition 
7P0 << ft and the TOD induced by the grating /3| is 
negligible, we can reduce Eqs. (1-2) to a single effective 
NLS equation. If we consider the evolution of the pulse 
in an NLFBG amplifier in the absence of gain saturation 
and for pulses with spectral bandwidth less then the am- 
plifier bandwidth, the pulse evolution can be described 
by an effective NLS equation with gain (NLSG) [7]: 
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where /3| and j g are the effective GVD and the ef- 
fective nonlinear parameter of the grating, g the gain 
coefficient and U is the Bloch-wave envelope. Since we 
know from previous analyses that the parabolic simi- 
lariton is solution to this equation, we could expect the 
pulse evolution in NLFBGs with gain to share common 
properties with classical similaritons. The term Bragg 
similariton refers specifically to the self-similar solution 
of this effective NLSG equation. 

In order to get a substantial increase of the effective 
GVD in an NLFBG operated in transmission, we need 
to propagate the pulse at a frequency very close to the 
photonic bandgap and/or use a grating with a very large 
index variation. In such cases, the effect of TOD cannot 
be neglected anymore and we have to evaluate how the 
BS is affected by higher-order dispersion. Fig. 1 shows 
the different dispersion coefficients for NLFBGs used in 
transmission. In these conditions, the GVD parameter 
(/?f) is many orders of magnitude higher than that of 
a classical fiber, but its TOD (/3f ) is also very high. It 
is now well known that the presence of TOD, or more 
precisely a large ratio, has a detrimental effect on 

the self-similar regime [8,9]; hense we have to minimize 
this ratio in order to generate a BS. We can also take note 
that the effect of TOD will be minimized for a narrow 
pulse spectrum. The fourth-order dispersion is also 
non-negligible and, as we will see, it has a significant 
effect on pulse propagation near the bandgap. 



, Group-Velocity Dispersion (GVD) , Third-Order Dispersion (TOD) 

i J , , , , , , , , in J , , , , , , , 




5 10 15 20 25 30 35 10 *0 5 10 15 20 25 30 35 



Fourth-Order Dispersion (FOD) A? /Pi ratio 
\ 2 r , , , , , , 3 n sn, . , , , 




Detuning [nm] Detuning [nm] 



Fig. 1. Dispersion coefficients for an NLFBG in trans- 
mission near the bandgap for three different values of 
index modulation: Sn = 10 -2 (blue), Sn = 5.0 x 10 -2 
(green) and Sn = 10 -3 (red). 

Following the notation of Bale et al. [8] we character- 
ized the effect of the higher-order dispersion coefficients 
on the propagated pulses in the NLFGB with the skew 
parameter which quantifies the asymmetry of the pulse 
temporal shape: 

S = J A 2 Z 3 d£ (4) 
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Fig. 2. Evolution of the skew parameter and parabolic 
misfit with detuning for an input gaussian pulse of 1.0 
(blue), 3.5 (green) and 6.0 nJ (red). 

where £ is a normalized time variable and A the nor- 
malized amplitude. The skew parameter has been com- 
puted for different values of detuning by propagating a 
5.0-ps gaussian pulse over a distance of 1 m with a clas- 
sical Split-Step Fourier (SSF) algorithm with the corre- 
sponding higher-order dispersion (P n <h) and the follow- 
ing NLFBG parameters: Sn = 1.0 x 10" 2 , g = 29 dB/m, 
7 = 54 km -1 W 1 . The misfit parameter between the 
pulse intensity profile \U\ 2 and the parabolic fit \Ufu\ 2 
has also been computed for the same conditions: 

ej = J (\U\ 2 - \U fu \ 2 fdt/ J \U\ A dt (5) 

Fig. 2 shows the skew and the misfit parameters (e 2 ) 
for various values of detuning and initial pulse energy. 
Interestingly, if we are assuming that these two param- 
eters are good indicators of the departure of the pulse 
from the classical similariton, this departure increases 
with detuning even if the ratio is decreasing. In 

the present case, the optimal detuning for the formation 
of a similariton-like pulse is between 6-8 nm. The gray 
region in the upper graph of Fig. 2 shows the condi- 
tion for which a shock- wave type of instability normally 
starts to develop (S > 0.06) in the absence of FOD. 
However, in the present case, the large value of FOD in- 
hibits the formation of instabilities, even in the case of 
a strongly distorted pulse. Fig. 3 shows some examples 
of classical similaritons distorted by TOD with differ- 
ent skew parameters with and without FOD. The pulse 
has been propagated over four values of distance with a 
classical SSF algorithm using the dispersion parameters 
/3f = 24.7 ps 2 /m and /?f = 4.5 ps 3 /m. We can clearly 
see the onset of wave breaking on the left side of the 
pulse propagating without FOD when, for the same con- 
ditions, the pulse propagating with /3| = 1.1 ps A /m is 
immune to this kind of instability. 

Considering the optimal conditions previously found, 
we simulated pulse propagation in an NLFBG with gain 
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NLFBG with gain in the normal dispersion regime. The 
stabilizing effet of FOD on the similariton propagation 
in the presence of TOD has been demonstrated. Finally, 
we have shown that pulses amplified in an NLFBG share 
common properties with classical similaritons generated 
in standard optical fibers with gain. This solution could 
open the possibility of parabolic amplification at wave- 
lengths usually confined into the anomalous dispersion 
regime such as Er-doped fiber amplifiers. 
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gineering Research Council of Canada. M. Laprise e-mail 
address is martin. laprise. l@ulaval.ca. 



Fig. 3. Effect of fourth-order dispersion on the wave 
breaking on an asymmetric pulse for different values of 
skew parameter: S « 0.05 (a), S « 0.06 (b), S « 0.07 (c) 
and S « 0.13 (d), with (blue) and without FOD (green). 



by solving numerically Eqs. (1-2) to all orders of dis- 
persion with the Optimized Split-Step Method (OSSM) 
described in [10]. The NLFBG used in these simulations 
has a length of 1 m, a constant gain coefficient g=29 
dB/m, a coupling constant ^=24.5 cm -1 and a detun- 
ing constant £=-60.5 cm -1 . The dispersion coefficients 
are given by: /3f = 83.4 ps 2 /m, /3| = 23.9 ps 3 /m and 
/3f = 9.8 ps 4 /m. The effective refractive index in the 
absence of the grating was equal to n=1.45 and the non- 
linear coefficient equal to 7 = 54 /cm _1 VK _1 . The input 
gaussian pulse has a temporal width of 5.0 ps and a peak 
power of 1.13 kW. In order to compare the results with 
propagation in a simple medium, the simulations have 
also been repeated with the classical SSF algorithm with 
the corresponding dispersion and nonlinear coefficients. 
Fig. 4 shows the temporal pulse characteristics on lin- 
ear and logarithmic scales from both methods. Since the 
OSSM is a temporal method and the group velocity in 
the NLFBG is changing with detuning, the three tem- 
poral profiles are identified by their propagation time. 
The profiles obtained with the OSSM and the classical 
SSF algorithm are very similar. We can suppose that 
the minor differences between the two curves originate 
from the cumulative effects of the higher-order disper- 
sion coefficients (/3 n >4) which were left out in the SSF 
simulations. More interestingly, the pulses share many 
properties with classical similaritons. The central part 
of the pulse evolves from a gaussian profile to a slightly 
skewed parabolic shape. When plotted on a logarithmic 
scale, the intensity profile presents abruptly decreasing 
wings typical of compact parabolic similaritons. The in- 
tensity profile also presents exponentially decaying low- 
intensity wings typical of the intermediate asymptotic 
regime of similaritons. Finally, the chirp in the high in- 
tensity region of the pulse exhibits an excellent linearity. 

In summary, we have numerically investigated the op- 
timal condition for the amplification of a pulse in an 
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Fig. 4. Temporal profile after 1635 (blue), 3272 (green) 
and 4906 ps (black) of propagation in the NLFBG using 
a classical SSF algorithm (a) and the OSSM (b). The 
gray curve in (a) corresponds to a classical similariton 
without TOD nor FOD. 
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